Multipole interaction between atoms and their photonic environment 
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Macroscopic field quantization is presented for a nondispersive photonic dielectric environment, 
both in the absence and presence of guest atoms. Starting with a minimal-coupling Lagrangian, 
a careful look at functional derivatives shows how to obtain Maxwell's equations before and after 
choosing a suitable gauge. A Hamiltonian is derived with a multipolar interaction between the guest 
atoms and the electromagnetic field. Canonical variables and fields are determined and in particular 
the field canonically conjugate to the vector potential is identified by functional differentiation as 
minus the full displacement field. An important result is that inside the dielectric a dipole couples to 
a field that is neither the (transverse) electric nor the macroscopic displacement field. The dielectric 
function is different from the bulk dielectric function at the position of the dipole, so that local-field 
effects must be taken into account. 

PACS numbers: 41.20.Jb, 42.50.-p 



I. INTRODUCTION 

Optical properties of atoms such as spontaneous- 
emission rates can be strongly influenced by their dielec- 
tric environment 0. It is well known that near a mirror, 
emission rates can be enhanced or diminished, depending 
on the distance to the mirror and the orientation of the 
atomic dipole moment Inside optical cavities, life- 

time effects are even stronger Q. In three-dimensional 
photonic crystals with a photonic band gap, spontaneous 
emission would even be fully inhibited at any position 
in the crystal for atomic transition frequencies within 
the band gap [5|. Not only single-atom properties, but 
also properties of two or several of atoms such as dipole- 
dipole interactions and superradiance will be influenced 
by the dielectric environment. Again, extreme changes 
compared to free space can occur for atoms positioned 
inside microcavities 0, 0] or photonic crystals [1, ■ 

The above medium-modified processes must be de- 
scribed by a quantum optical theory of dielectrics. In 
this paper such a theory will be given of guest atoms in- 
teracting with a photonic dielectric environment that is 
characterized by a given spatially varying and real dielec- 
tric function £(r). The guest atoms by definition are the 
atoms which are not included in the dielectric function. 
These guests will be described microscopically, whereas 
the dielectric is described macroscopically in terms of the 
dielectric function only. 

Dielectric mirrors and photonic crystals are usually 
described by a frequency-independent spatially vary- 
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ing refractive-index. Optical components such as glass 
plates, lenses and optical cavities are some more exam- 
ples. In many cases, the refractive index can be con- 
sidered as piecewise constant, but not always: in the 
so-called graded-index optical fibres the refractive index 
in the core varies parabolically with the radius I'^ 
this paper, the relative dielectric function £(r) is left un- 
specified (but assume it to be piecewise continuously dif- 
ferentiable) so that the theory describes both piecewise 
constant and continuously varying dielectric functions. 

The quantization of the electromagnetic field in free 
space can be found in man y te xtbooks on quantum op- 
tics [Hill El El. In mill El El the more general 
problem is addressed how to quantize the electromag- 
netic field in a dielectric described by a given real dielec- 
tric function e(r) that depends on position. The term 
"macroscopic quantization" has been coined for this pro- 
cedure E3 ■ The special case of infinite photonic crystals 
is treated in E3- 

Guest atoms can be described theoretically essentially 
in two ways. The simplest way is to treat them as 
known probes of the electromagnetic field in the medium. 
In that case, the electromagnetic fields are found from 
Maxwell's equations in the absence of the guest atoms; 
the guests are introduced as atoms with given proper- 
ties, such as transition frequencies and dipole moments. 
The atoms are assumed to couple to the fields that were 
found in their absence. The second and more funda- 
mental way to introduce guest atoms into the theory is 
to start with Maxwell's equations that also contain as 
sources the charges that make up the guest atoms. In 
this second approach, a Hamiltonian for the combined 
system of charges and fields should be found that leads to 
Maxwell's equations, both inside and outside the atoms. 

Even for an atom in free space, the difference between 
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these two approaches has led to debates how to interpret 
the field to which a dipole couples, either to the trans- 
verse part of the electric field or to the displacement field 
113, lH m mil. The latter couphng is the correct 
(more fundamental) interpretation 13], but for most ob- 
servables there are no numerical punishments when in- 
terpreting the field wrongly. However, inside a dielectric, 
there would be a considerable difference between an atom 
coupling to the electric field or to the displacement field. 
In this paper it will be shown by using the second, more 
fundamental way of introducing guest atoms, that nei- 
ther dipole coupling is correct in a dielectric. Moreover, 
the need to consider local-field effects will arise in a na- 
tural way. 

Of course, one does not tell the whole truth about a 
dielectric when describing it by a real, nondispersive and 
spatially varying dielectric function £(r). In a macro- 
scopic description one forgets details of the microscopic 
constituents of the dielectric. Also, material dispersion 
and absorption of light (transitions to nonradiative states 
in the dielectrics) are neglected. It is well known that 
the dielectric function is a response function that should 
be a dispersive and complex function of frequency, so 
as to satisfy the Kramers-Kronig relations. Certain sum 
rules ii2&] for modified spontaneous emission rates when 
averaged over all frequencies will therefore not hold in 
the present formalism that violates these relations. The 
question is whether dispersion and absorption are im- 
portant in a particular experiment that one has in mind. 
Often optical experiments are only interesting in a fre- 
quency range where absorption is indeed negligible, for 
example when measuring light emitted by excited atoms 
inside a photonic crystal. In such cases, it is common 
practice to neglect material dispersion and absorption in 
the theoretical description as well 0, . 

Quantum optical descriptions exist where dielectric 
functions do satisfy the Kramers-Kronig relations, both 
for homogeneous l23j and inhomogencous dielectrics 
|29, 30, 31, 32). Usually, the guest atoms are introduced 
into these theories in the simplest of the two ways de- 
scribed above, as probe atoms with known properties in 
a two- or three-level description. It would be interesting 
to introduce guest atoms in theories of inhomogencous 
Kramers-Kronig dielectrics in the more fundamental way, 
starting with Maxwell's equations with the charges of the 
guest atoms as sources, but this will not be done here. 
Nor will we look at more microscopic descriptions of the 
dielectric 1^3, 0, 113 where dispersion shows up natu- 
rally. It would be very challenging to derive optical pre- 
dictions from a microscopic description of light and the 
matter that builds up a complex dielectric such as a pho- 
tonic crystal. 

The goal of this paper is to derive a Hamiltonian 
with multipolar interaction between the guest atoms and 
the electromagnetic field inside the inhomogencous and 
nondispersive dielectric. The dipole Hamiltonian can 
then be found as an approximation. First, field quan- 
tization of a dielectric without guest atoms is described 



in Sec. [HJ Atoms are introduced into the dielectric in 
Sec. mil The starting point will be Maxwell's equations 
and a minimal-coupling Lagrangian that produces these 
equations. Special attention is paid to check whether 
Maxwell's equations still hold after choosing a gauge. 
This requires an interesting analysis of functional differ- 
entiation after choosing a gauge, presented in Sec. IIVI 
The minimal-coupling Lagrangian is then transformed in 
Sec. to its multipolar form with use of the Power- 
Zienau-WooUey transformation that is well-known for 
free space careful analysis of functional differ- 

entiation allows us to use a transformation that is simpler 
and more like the free-space case than presented in re- 
lated work 0, . After the transformation, canonical 
variables and fields are determined. In particular, the 
important question which field in the dielectric is canon- 
ically conjugate to the vector potential can be answered 
more easily than in fvA IT^ and the answer will be dif- 
ferent than presented in |17l |. The multipolar Hamilto- 
nian and its dipole approximation are given in second- 
quantization notation. Results are compared to the free- 
space case. The free-space dipole-coupling controversy 
and confusion is reviewed in Sec. IVII and the results of 
this paper are discussed in that perspective. The quan- 
tum optics of dielectrics which also have inhomogencous 
magnetic properties are briefiy discussed in Sec. lVIII be- 
fore concluding in Sec. IVIIII 

II. INHOMOGENEOUS DIELECTRIC 
WITHOUT GUEST ATOMS 

A. Classical Lagrangian and Hamiltonian 

In this section the quantization of the electromagnetic 
field in inhomogencous dielectrics 0, is reviewed. 
The emphasis will be on concepts and results that will 
be employed in the following sections, when guest atoms 
are introduced in the dielectric. 

In Sl-units, the source-free Maxwell equations in mat- 
ter are 

V-B==0, VxE + B = 0, (la) 

V-D==0, VxH-D==0. (lb) 

Here and in the following the dot denotes a partial time 
derivative. The fields E, B, D and H are the electric field 
and the magnetic induction, the displacement field and 
the magnetic field vector, respectively. For nonmagnetic 
inhomogencous dielectrics, the constitutive relations are 
simply B — hqH and D = eoe{r)'E. (The generalization 
of the present theory to inhomogencous magnetic mate- 
rials will be discussed in section IVHI ') As for free space, 
the electric and magnetic fields can be expressed in terms 
of a vector potential A and a scalar potential <I>: 

E=-V$-A, B = VxA. (2) 

There is gauge freedom in choosing pairs (A, $) that lead 
to the same electric and magnetic fields. Now choose 
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the generahzed Coulomb gauge which is defined by the 
requirement that the vector potential satisfies 



V • [e(r)A(r) ] 0. 



(3) 



The vector potential or any field satisfying this condi- 
tion, is called "generalized transverse" , because it satis- 
fies a generalized version of the Coulomb gauge condi- 
tion V • A = in free space. In the generalized Coulomb 
gauge, the vector potential must satisfy the wave equa- 
tion 



e(r) .. 

V X V X A + At^A = 0, 



(4) 



in order to be consistent with the Maxwell equation l|lbp . 
The scalar potential can be chosen identically zero ($ = 
0) in the generalized Coulomb gauge. 

Since the goal is to find a quantum optical Hamiltonian 
in the end, one should start with a Lagrangian formalism 
for the classical Maxwell fields. From the Lagrangian 
the canonical fields and their conjugates can be identified 
that will become pairs of non-commuting field operators 
in a later stage. The principle of least action states that 
fields (and particle variables, when present) minimize the 
action the action is defined as the time- integrated 
Lagrangian between some initial and final times. The 
requirement that small variations in the fields do not 
change the action leads to the Euler-Lagrange equations 
for the canonical fields, in our case for the vector poten- 
tial [mil: 



6L d SL _^ 
SA dtJX ^ 



(5) 



and 



Here, functional derivatives are denoted with "J" 
more will be said about them later. 

A Lagrangian Lq for the electromagnetic field in an 
inhomogeneous medium is 



Lo 



dr Cq 



dr 



MoHVx A)2 



(6) 

The Lagrangian is the spatial integral of the Lagrangian 
density Co over a large volume V that will eventually 
be sent to infinity. The subscript "0" is used to denote 
the absence of guest atoms in the dielectric. The vector 
potential is a canonical field variable and its canonically 
conjugate field can be found as a functional derivative of 
the Lagrangian density 
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SA 



eoe(r)A = -D. 



(7) 



(The functional derivative is used somewhat naively here, 
but the answer is correct, as a more detailed analysis in 
section Hvl will show.) In other words, the field canoni- 
cally conjugate to the vector potential equals minus the 
displacement field, which is a transverse field. Proceeding 
as for free space [13, one finds that the Euler-Lagrange 
equation of motion for the vector potential leads to the 



wave equation Q for the vector potential in the medium. 
The Hamiltonian is 



Ho= / dr (n.A-£o) = ^ / dr 



n2 



(V X A)' 



eoe{r) 



Mo 



(8) 

This is the Hamiltonian for the classical electromag- 
netic field in an inhomogeneous dielectric, without guest 
atoms. 



B. Complete sets and quantum Hamiltonian 

For a quantum optical description of the dielectric, 
the electromagnetic fields can best be expanded in terms 
of harmonic solutions of the wave equation With 
each of these "true modes" one can associate independent 
canonical variables, for which commutation rules can be 
given. The set of true modes is not unique. This freedom 
will be used below to choose a particularly convenient 
set. For example, in vacuum the true modes are trans- 
verse plane waves. For the plane waves, one can choose 
linear combinations of cosine and sine solutions cos(k- r) 
and sin(k-r). The complex exponential exp(ik-r) is only 
one such linear combination. 

It will now be argued why the field can be expanded in 
terms of a set of real mode functions, and the orthonor- 
mality relations of the modes will be derived. The choice 
of real modes will simplify the quantization procedure. 
The reason not to start with complex mode functions is 
that the latter has associated generalized coordinates and 
momenta that are not Hermitian. This makes interme- 
diate results more complicated 0, 0, [s^ . Real mode 
functions were also used in [T5l |. without the motivation 
given here. 

Why does a complete set of real mode functions exist? 
As before, the electromagnetic fields are assumed to live 
in the large volume V. Let Q be an abstract operator in 
Hilbert space which has a local representation (r|Q|r') = 
(5(r — r')Q(r) in configuration space: 



Q{r) 



= V X V X 



(9) 



The operator Q is Hermitian under the normal inner 
product. Eigenvalues of Q are [uox/cf and the uj\ will 

be called eigenfrequencies^ All eigenfunctions gA of Q 

have the property V • [\/£(r)g(r)] = 0. The label A is 
understood to count both continua and discrete sets of 
solutions. The subspace of functions in Hilbert space 
with the same transversality property is spanned by the 
eigenfunctions of Q. Now let C be the operator which 
is also local in configuration space, where its action is 
to take the complex conjugate. The dielectric function 
in this context must also be viewed as an abstract op- 
erator £ with local representation in configuration space: 
(r|£|r') = 5{y — r')£(r). The representations of Q and 
C commute in configuration space, because £(r) is real 
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for all positions r in y. Then Q and C commute in any 
representation. From the fact that Q and C commute it 
follows that an orthonormal basis of real eigenfunctions 
{gA} of Q can be chosen to span the subspace (with com- 
plex coefficients). 

The above analysis shows that the vector potential A 
can be expanded in terms of a complete set of real vector 
mode functions {hA(r)} = {■y/e(r)gA(r)}, which are the 
harmonic solutions of the wave equation (0J: 

V X V X hA(r) - ^^hA(r) = 0. (10) 

These mode functions satisfy the same generalized 
transversality condition (PJ as the vector potential. As 
is clear from Eq. (|10|l , unlike the gA the functions \v\ do 
not satisfy a Hermitian eigenvalue equation. (To each 
typ e of modes corresponds a different density of states 
[STj l . ) From the orthonormality of the gA it follows that 
the functions hA(r) satisfy the generalized orthonormal- 
ity condition ,1^] 

ydre(r)hl(r).hA'(r) = 5AA'. (11) 

The complex-conjugation symbol * was written for future 
reference, since of course \\\ is real. The equation ifTT)) 
will be called a generalized inner product of the modes 
Ha and . The spatial integral in Hll|l runs over the 
volume y, so that the mode functions scale as (y)"^/^. 

In a scattering situation, where e(r) is a space-filling 
dielectric function plus a local modification within a scat- 
tering volume , the contribution of the scattering vol- 
ume to the integral (|ll|l scales as 14 /l^- This fraction 
becomes of measure zero when the quantization volume 
V is sent to infinity. As an example, mode functions of an 
infinite photonic crystal with a single point defect (extra 
or missing dielectric material) have the same orthonor- 
mality relations as the mode functions in the absence of 
the defect. 

The functions \v\ are complete in the sense that they 
form a basis for generalized transverse functions (such as 
the vector potential) that satisfy the wave equation Q). 
In other words, a generalized transverse delta function 
(a distribution) can be defined in terms of the functions 
Ha: 

5j(r,r')^5]hA(r)hA(r')e(r'). (12) 

A 

For £(r) = 1, this expression reduces Jo the (real) free- 
space transverse delta function (see [l^, p. 53) 

<5T(r)^^5(r)l-^(l-3f0r), (13) 

where I is the unit tensor in three dimensions and r is the 
unit vector r/|r|. Evidently, the generalized transverse 
delta function (|12|l is real because the mode functions 
are real. From its definition (|12|l and the normalization 



condition (|ll(l of the modes, it follows that is idem- 
potent: 

1 dn 5j(r,ri)-^J(ri,r') =<5j(r,r'), (14) 

In other words, (5 J is a projector into the subspace of gen- 
eralized transverse functions. The generalized transverse 
delta function is not symmetric in its arguments, because 
it is transverse in its second and generalized transverse 
in its first variable: 

Vr • [£(r)jj(r,r')] = = Vr^ • [^J(r,r')] . (15) 

The vector potential and its canonically conjugate field 
have normal-mode expansions 

A(r,0 = l/(V^)^<7A(i)hA(r), (16a) 

A 

n(r,i) = V^^pA(i)e(r)hA(r), (16b) 

A 

with generalized coordinates <i\(t) and momenta p\{t). 
At this point the choice of real mode functions pays off, 
because the associated generalized coordinates and mo- 
menta are only real when the modes are real; only real 
coordinates and momenta will become Hermitian oper- 
ators in a quantum description. If the above two ex- 
pansions are substituted in the Hamiltonian ((SJ and the 
normalization condition (|ll|l is used, then it follows that 
H — \ X]a(Pa +'^a9a)- The Hamiltonian turns out to be 
a simple sum over the true modes of the inhomogeneous 
dielectric, where the energy of each mode corresponds to 
a one-dimensional harmonic oscillator with position qa, 
momentum p\ and frequency uj\ . 

Now comes the quantization step. The independent 
canonical pairs satisfy the standard equal-time commu- 
tation relations [q\{t)TP\i{t)] = ih5\y. With the normal 
mode expansions ((TCajl and lfT6B)l . the commutation rela- 
tion for the vector potential and its canonically conjugate 
field can be found immediately: 

[A(r,i),n(r',t)] = *ft^hA(r)hA(r')e(r') =*?M5j(r,r'). 

(17) 

The commutator is a dyadic quantity. It turns out to be 
proportional to the generalized transverse delta function. 
Annihilation operators are introduced as a\ = 
{2h)q\ + iyj\/{2fujj\)p\ and creation operators as 
their Hermitian conjugates. They have standard com- 
mutation relations [a\{t), a\,{t)\ = (5a a' and all other in- 
equivalent commutators are zero. The Hamiltonian be- 
comes the sum over contributions hijjx{a\a\ -\- 1/2) of 
individual modes. Thus the concept of a photon as the 
elementary excitation a^|0) of a mode is as useful for in- 
homogeneous dielectrics as it is for free space. Number 
states, coherent and squeezed states etcetera can be de- 
fined analogously. The only difference for inhomogeneous 
dielectrics is that their true modes are not plane waves. 
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The vector potential operator and its canonically con- 
jugate field operator can be expressed in terms of creation 
and annihilation operators as 



A(r,i) = J2 



hA(r)+H.c.], (18a) 



U{r,t) = -zeo£(r)^./|^[aA(t)hA(r)-H.<48b) 



2eo 



Here, "H.c." denotes the Hermitian conjugate. The 
forms of the electric and magnetic fields as quantum 
mechanical operators as well as their commutation rela- 
tions immediately follow from ^ and the above equation 
(|18a|l . The time-dependence of the operators is simply 
harmonic, for example a\(t) = a(0) exp^—iujxt) . 

In practice, it can be convenient to use a set of complex 
true mode functions {f;^} instead of the real mode func- 
tions {h^}. Since the complex mode functions should 
also satisfy the wave equation (|10() and the general- 
ized orthonormality condition Hll|l . the two sets of mode 
functions are related through a unitary transformation 

= U^x hx that only relates mode functions with 
identical eigenfrequencies; U is a unitary matrix. Note 
that because of this unitary relation, the generalized 
transverse delta function H12|) can alternatively be ex- 
pressed in terms of the complex mode functions f^. Its 
effect is the substitution of one by ix and the other 
hA by i*^: <5j(r,r') is also equal to fA(r)f*(r')e(r'). 
However, after this substitution it is no longer obvious 
that Sj{r,r') is real- valued. 

The field operators can also be expanded in terms of 
the complex mode functions fx as 



U{r,t) = -ieae{r)Y^ 



hhjx 
2en 



ai°)(i)fA(r)-H.(=lf)b) 



where the new annihilation operator a^'^^ associated with 
the complex mode f^^ is defined in terms of the "old" 

operators as afi!^ = ax, and o'^^'^ is its Hermi- 



tian conjugate. The commutation relations of a^l^^ and 
are again the standard relations, because U is a uni- 



JO)t 



tary transformation. To distinguish a^^^ from operators 
to be defined later, the superscript (0) has been added, 
signifying that no guest atoms are present. The time 
dependence of the operators is again harmonic. 

This completes the quantization of the electromagnetic 
field in an inhomogeneous dielectric without guest atoms. 
The reason to start the quantization procedure with real 
mode functions was that the associated generalized co- 
ordinates and momenta are real quantities that become 
Hermitian operators in quantum mechanics. It is possible 
to start with complex mode functions instead and to pro- 
ceed with the non-Hermitian operators ESimil, but 



it makes intermediate results unnecessarily more compli- 
cated. The unitary relations between complex and real 
mode functions and between their respective annihilation 
operators are purely formal, unless both sets of mode 
functions are given explicitly. The above quantization 
procedure only relies on the mere existence (rather than 
on an explicit construction) of these unitary mappings. 



III. INHOMOGENEOUS DIELECTRIC WITH 
GUEST ATOMS 



In the previous section it was described how to quan- 
tize the electromagnetic field in an inhomogeneous di- 
electric. Now inside the inhomogeneous dielectric guest 
atoms are introduced. Their optical response is not in- 
cluded in the dielectric function e(r) of the medium. The 
goal in the following sections is to find the quantum op- 
tical description of the combined system, with a mul- 
tipole interaction between the electromagnetic field and 
the guest atoms. There are at least two reasons why 
the multipolar Hamiltonian is to be preferred. In the 
first place, it is more convenient when only approximate 
calculations can be done which in the minimal-coupling 
formalism would give gauge-dependent results 11]; sec- 
ondly, atoms are much smaller than optical wavelengths 
and in the multipole-formalism this can be exploited well. 
Actually, atoms are so much smaller than optical wave- 
lengths that often "atoms" are identified with "dipoles" . 

The starting point is the minimal-coupling Lagrangian 
that produces the Maxwell equations and the equations 
of motion for the charges that make up the guest atoms. 
The minimal-coupling Lagrangian can be used to find 
a minimal-coupling Hamiltonian and this procedure can 
be found in [l^ [l^ . A clear exposition is also given in 
|38l |. Here the Lagrangian will first be transformed to 
the multipolar form before constructing a Hamiltonian. 
The latter procedure was followed also in 0, 0| . The 
present work is different in some essential aspects that 
will be stressed where appropriate. 



A. Choice of suitable Lagrangian 

Guest atoms inside an inhomogeneous dielectric can be 
described by a charge density Cg and a current densit y 
which show up as sources in Maxwell's equations pTl[l7l |: 



V-B = 0, 

V X E + B = 0, 

£oV • [e(r)E(r)] = 

^0 X B - £o£(r)E = Jg 



(20a) 
(20b) 
(20c) 
(20d) 



Here, Ug is the charge density and Jg the current density 
produced by the guest atoms alone, as stressed by the 
subscript "g" ; the dielectric is completely described by 
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the dielectric function eoe{r) and the magnetic perme- 
abihty /ig. Whatever Lagrangians and Hamiltonians are 
introduced for the inhomogeneous dielectric plus guest 
atoms, they must lead to these four Maxwell equations. 
Moreover, the electrons with charges — e and masses rric 
should respond to electric and magnetic fields as given in 
the equation of motion 



[E(r„j) + Tmj X B(r™j)] 



(21) 



We assume that there are no free charges. All electrons 
(labelled j) are bound to atomic nuclei (label m) to form 
neutral guest atoms. Then dg and Jg are given by |39| 



(r,i) 



Jg(i-,i) = -eY^ 



E/ ^"'■J 



j 



(22b) 



Here, is the nuclear charge of atom m. The guest 
atoms are assumed to have fixed positions, their nuclei 
are their centers of mass and are stationary at posi- 
tions Rm. From these explicit forms of Ug and Jg fol- 
lows the equation of continuity or current conservation, 
V-Jg-t-cTg = 0, which can also be found from the Maxwell 
equations (|20cl) and H20dp . 

Again, the electric and magnetic fields can be defined 
through Eq. Q in terms of a vector potential A and a 
scalar potential $. Then the two homogeneous Maxwell 
equations H2Ua|) and Ij20bp are automatically satisfied. 
The other two Maxwell equations should follow from the 
Euler-Lagrange equations [Eq. for the scalar and the 
vector potential, respectively. The minimal-coupling La- 
grangian is 



dr Lr, 



(23) 



Here the Lagrangian density £inin describes the electro- 
magnetic field energy and its minimal-coupling interac- 
tion with the guest atoms: 



-Cmin = ^eoe(i") 



1 



A + V$ -- — (Vx A)^+Jg-A-crg$. 
J 2/^0 

(24) 

Indeed, Maxwell's third and fourth equations can be 
found from the Euler-Lagrange equations for the scalar 
and the vector potentials, respectively. Moreover, the 
Euler-Lagrange equations for the canonical variables Vmj 
give the equations of motion (|21|l for the charged parti- 
cles. Note that the Lagrangian leads to these equations 
of motion, before choosing a gauge to fix A and <i> with: 
the equations of motion are gauge-independent results 
that should not depend on the choice of gauge. 



B. Fixing the gauge 

The electric and magnetic fields are defined in terms 
of a scalar and a vector potential. But there is gauge 
freedom, which means that the scalar and vector poten- 
tials are not uniquely defined by the requirement that 
measurable electric and magnetic fields satisfy Maxwell's 
equations. We need to choose a gauge in order to find in 
the end a quantum mechanical description of light inter- 
acting with the guest atoms. As in the situation without 
guest atoms in section^ the generalized Coulomb gauge 
is chosen so that the vector potential satisfies Eq. . In 
this section it will be checked whether the equations of 
motion for the scalar and vector potentials still lead to 
the third and fourth Maxwell equations after choosing 
the gauge. This must be the case, because the choice of 
gauge should not change the physical predictions of the 
theory. Still, the check was not performed in [TtLIi^ and, 
as we shall see, it will be very useful to do so here. 

The gauge affects the interaction term 



La$ = £o y dr e(r)A • V*, 



(25) 



of the Lagrangian (|23|l . The term becomes identically 
zero, because in the generalized Coulomb gauge it has 
become an inner product of a transverse and a longitu- 
dinal function. The remaining terms in the Lagrangian 
involving the scalar potential lead to an Euler-Lagrange 
equation that is the generalized Poisson equation for the 
scalar potential in the Coulomb gauge: 



eoV • [e(r)V$(r)] = -a^. 



(26) 



Clearly, the scalar potential can not be chosen identically 
zero as in the situation without guest atoms. The gauge- 
fixing condition (jS)) for the vector potential, the equation 
(|26|l for the scalar potential, together with the definition 
of the electric field Q in terms of the two potentials, still 
lead to the third Maxwell equation (|20c|) . 

As in free space, the scalar potential is a function of 
the positions of the charges that make up the guest atoms 
|T3l |. In other words, one can first solve the coupled equa- 
tions of motion for the vector potential and the charges, 
and from the charge distribution (Jg{t) thus found, the 
scalar potential $(t) can be found as the solution of 
Eq. H2t)|l . Therefore, the scalar potential is not an inde- 
pendent canonical field. The Lagrangian H23|l can then 
be simplified as 



i'min = XI " 



dr C'r 



(27) 



where Eqs. 1)25(1 and 126(1 were used. The Coulomb in- 
teraction Vc = (so/2) / dr e(r)(V<I')^ is a function of 
the guest atoms alone; the Lagrangian density in 1(27(1 
becomes 



1 



1 

2/io 



(V X A)2 + Jg . A. 



(28) 
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Which equation do we find for the vector potential after 
choosing the generahzed Coulomb gauge? Let us begin 
at the other end: in order to be consistent with the fourth 
Maxwell equation Ij20d|l . the vector potential should sat- 
isfy 

V X V X A + eo£(r)A = Jg - eoe(r)V$. (29) 

It is not obvious how the source term — £o£(r)V<I> can 
appear at the right-hand side of this equation by func- 
tional differentiation of the Lagrangian with respect to 
the vector potential. Before choosing the gauge, this 
source term originated from the interaction term Xa* 
[Eq. (|25|l ]. which is zero after choosing the gauge. After 
choosing the generalized Coulomb gauge, a more careful 
analysis is needed in order to find the fourth Maxwell 
equation. 

IV. FUNCTIONAL DIFFERENTIATION AFTER 
CHOOSING THE GAUGE 

After choosing the gauge, the vector potential is gen- 
eralized transverse. The Euler-Lagrange equation for the 
vector potential is therefore an equation of motion of a 
constrained system, where the constraint is the gauge 
condition ||2Jl . One could try and solve this problem using 
the method of Lagrange multipliers 0| , but this is not 
the route that will be pursued here. Instead, the appro- 
priate mathematical definition and computation of func- 
tional derivatives after choosing a gauge will be studied 
in the following subsections II V Al and II V 51 respectively. 
The results will be applied to our physical problem in 
subsection IIVCI 

A. Two definitions of functional derivatives 

Let us generalize the problem somewhat by considering 
a functional F = J dr !F with a functional density T that 
depends on the three-dimensional vector fields X and Y. 
Assume also that at some stage the generalized Coulomb 
gauge will be chosen for the field Y. This gauge is defined 
by the requirement that Vr • [e(r)Y(r)] equals zero. 

Before choosing the gauge, the functional derivative of 
the functional F with respect to the vector field Y is 
defined as 

SF fdr'{.F[Y(r') + 7^(r-r')l]--^[Y(r')]} 

- — -— = hm 

SY(r) 7^0 7 

(30) 

(The X-dependence of was dropped for brevity.) The 
functional derivative of F with respect to Y describes the 
relative changes of F when small variations proportional 
to 6{r — r')l are added to the vector function Y. Here, 
(5(r — r') is the Dirac delta function in three dimensions; 
as before, I is the unit tensor. It turns out that the right- 
hand side of Eq. H3U|I can be computed as the partial 
derivative of the functional density J- with respect to Y. 



While doing this, the T can simply be considered as a 
function and Y as one of its variables. Before choosing 
the generalized Coulomb gauge, Eq. (pHHl correctly defines 
the functional derivative of F with respect to Y. 

Now suppose for the moment that F is defined as 
F = / dr X • Y. Suppose also that the field X is the 
product of e(r) with some longitudinal vector field. Then 
F becomes identically zero in the generalized Coulomb 
gauge, because it is the inner product of a transverse and 
a longitudinal vector field. Still, the functional derivative 
Eq. (|30|) of F would give a nonzero answer. This can only 
mean that Eq. I|3U|) does not define the functional deriva- 
tive with respect to generalized transverse functions cor- 
rectly. The reason is that the function space in which 
the field Y lives has become smaller by choosing the 
gauge: it now lives in the subspace of functions which are 
generalized transverse. This also means that functional 
variations of Y should stay inside this subspace. In the 
functional derivative (|30|l . variations in the whole func- 
tion space are allowed and clearly Vr' ■ [e{r')S{r — r')l ] 
is nonzero. 

In general, with every set of constraints a new func- 
tional derivative can be associated. Here, only the gauge- 
constraint will be considered that functions be general- 
ized transverse. Functional differentiation with respect 
to g eneralized transverse functions can be defined as (see 

m, p- 20) 

SF ^ ^.^ / dr' myjr') + 7 Sjjr', r)] - .F[Y(rO]} 
(5Yj(r) 7^0 7 

(31) 

with the generalized transverse delta function 5j as de- 
fined in Eq. H12I) . In this new functional derivative, the 
functional variations do stay inside the generalized trans- 
verse subspace, since Sj is the projector into the subspace 
and Vr' • [e(r')<5j(r',r)] = 0. The derivative will be 
called the "constrained functional derivative" in the fol- 
lowing. 



B. Simple rules to compute constrained functional 
derivatives 

Now the goal is to find simple rules to compute the con- 
strained functional derivative H31|) with respect to gen- 
eralized transverse functions, just like the normal func- 
tional derivative (|30|l can simply be calculated as a par- 
tial derivative. With that goal in mind, first some prop- 
erties of generalized transverse functions are derived. 

With every transverse function X"^(r) a generalized 
transverse function [X'^(r)/e(r)] can be associated. From 
section III Bl we know that the latter function has an ex- 
pansion in terms of generalized transverse modes hA(r), 
so that X"'"(r) can be expanded in terms of £(r)hA(r). 
This simple fact, in combination with Eqs. and 112|l . 
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leads to the following projection properties of dj: 

J dr'xT(r')-Jj(r',r) = X.^{r), (32a) 

e{r) J dr'(5j(r,r')-XT(r')/e(r') = XT(r), (32b) 

J dr'<5j(r,r')-XL(r') = 0, (32c) 

J dr'e(r')XL(r')-C(r',r) = 0, (32d) 

for any transverse function X""" (zero divergence) and lon- 
gitudinal function X'" (zero curl). 

With the use of Eqs. (|32a|l and Ij32d|l the functional 
derivative H31|l is simple in the following two important 
cases: 

I dr' XT(r') • Y(r') = X^{r), (33a) 
-^^^ I dr' e(r')XL(r') • Y(r') = 0, (33b) 

where X"^ and X'^ are arbitrary transverse and longitudi- 
nal functions, respectively. The second case (|33b|l makes 
clear that the constrained functional derivative of in- 
ner products of transverse and longitudinal fields indeed 
gives zero; the first case H33a|l shows that the partial- 
derivative-of-jT computation rule still gives the correct 
answers for inner products of Yj with transverse func- 
tions. 

How can the constrained functional derivative be cal- 
culated in the more general situation 



JYj(r) 



J dr' X(r') . Y(r'), 



(34) 



where X is a general vector function? It will now be 
shown that any vector field X can be decomposed such 
that the only two rules of computation needed are the 
simple cases (|33a|) and (j33bll . 

Given the vector field X, construct the scalar field 
(Tx = — V-X. Now find the potential X, given the "charge 
distribution" crx and the dielectric function eoe{r), from 
the following generalized Poisson equation: 



eoV- [£(r)Vx(r)] - -ax(r). 



(35) 



This is a well-known problem in electrostatics. There is 
a unique solution for x of this inhomogeneous problem, 
given the charge distribution and the boundary condition 
that the potential be zero at infinity. With the potential 
X thus found, define two vector fields Xi and X2 as 

Xi = X(r)-eoe(r)Vx(r), (36a) 
X2 = £o£(r)Vx(r), (36b) 

so that evidently X = Xi -I- X2. The vector field Xi is 
transverse by construction of the potential x; the field 



[X2/e(r)] is of course longitudinal. In summary, the fol- 
lowing theorem was proven: An arbitrary vector field can 
be uniquely decomposed into a part which after division 
by e(r) is longitudinal, and a transverse part. This the- 
orem is useful for evaluating the constrained functional 
derivative, because it leads to 



6Yj{r) 



J dr'X(r')-Y(r') ==X(r)-£o£(r)Vx(r 



(37) 



where the unique decomposition of X was used and the 
simple derivatives Eqs. H33all and Ij33b|l were applied to 
Xi and X2, respectively. The problem of computing a 
functional derivative with respect to a generalized trans- 
verse function has thus been reduced to a problem in elec- 
trostatics. Note that the constrained functional deriva- 
tive (|37|l produces a field that is always transverse. That 
transverse field is equal to the transverse part of X if X 
itself is transverse or if £(r) = 1. 

The unique decomposition (|36a|l and Ij36b|) of vector 
fields is a generalization of the Helmholtz theorem [4ll |. 
but the name "generalized Helmholtz theorem" was al- 
ready given to a slightly different statement , namely 
Every vector field Zi can be uniquely decomposed as the 
sum of a generalized transverse vector field Zi and a lon- 
gitudinal field 7i2- (The proof of this theorem in begs 
the question whether the part of the decomposition that 
is called longitudinal indeed has zero curl, but one can 
show that this is the case.) As a corollary of the decom- 
position (|36a|) and Ij36bp . a new and short proof can be 
given of the generalized Helmholtz theorem. The proof 
is simple: given Z, define X — e{r)Z. Then apply the 
previous unique decomposition to X, as in Eqs. I|36a|) and 
(|36H|l . Define Zi = Xi/£(r) and Z2 = X2/e(r). Then 
it follows that Z = Zi -I- Z2, where Zi is a generalized 
transverse and Z2 is a longitudinal field. This completes 
the proof. 



C. Functional derivatives of the minimal-coupling 
Lagrangian 

The definition of the constrained functional derivative 
and its computation rules can now be applied to our case 
of interest, where the functional is the Lagrangian i'min, 
Eq. H27II . and where the generalized Coulomb gauge ap- 
plies to the vector potential A. 

Before choosing the gauge, the "ordinary" functional 
derivative (|30|) of the Lagrangian of Eq. H24|l with re- 
spect to A leads to the fourth Maxwell equation Ij20d|l . 
as it should. After choosing the gauge, the interaction 
term La* [Eq. (|25|l ] in the Lagrangian becomes identi- 
cally zero. Its functional derivative with respect to A 
should also be zero. This is indeed the case, because 
the constrained derivative is the correct one to use 
rather than the ordinary functional derivative l|30|) after 
choosing the gauge. Note that the constrained functional 
derivative must also be used for free space after choosing 
the Coulomb gauge, with 6j equal to 6"^ , Eq. H13|l . For 
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free space the machinery of functional derivatives usually 
is not introduced and the derivative is taken implicitly, 
for example in [3 (p. 289). 

The naive calculation of the canonical field in Eq. {Tj) 
of section Hll can now be justified: 

I dr' e{v')k\v') = 2 £(r)A(r). (38) 
5Ai (r) J 

One can find this result by realizing that the functional 
on the left-hand side is a special case of Eq. (|33a|l with 
the fields X and Y equal to e(r)A and A, respectively. 

In the special case that the vector field X is the current 
density Jg produced by the guest atoms in the dielectric, 
current conservation implies that the scalar field aj^ as 
constructed in section IjIV B|) equals the time-derivative 
of the physical charge density ag . By the uniqueness of 
the solution of the generalized Poisson equation, the po- 
tential X must then be identified with the time-derivative 
of the physical scalar potential <i>. Therefore, the con- 
strained functional derivative of J dr Jg • A can now be 
computed as 

J dr' Jg(r') . A(r') = Jg(r) - £oe(r) V4>(r). 

(39) 

With this result, the Euler-Lagrange equation for the vec- 
tor potential from the Lagrangian (|27() precisely becomes 
the equation H29|) for the vector potential that we were 
looking for. Only by the careful computation of func- 
tional derivatives as presented in Sees. IIV Al and IIV Bl 
can one prove that the fourth Maxwell equation Ij20d|) 
holds also after choosing the generalized Coulomb gauge. 
Interestingly, before choosing the gauge, the source term 
— eoe(r)V$ in (|29|l came from the La<s> interaction term 
(|25|l in the Lagrangian. After choosing the gauge, how- 
ever, the source term is produced by the constrained 
functional derivative of the minimal-coupling interaction 
term / dr Jg • A. 

The left-hand side of the wave equation (^5)1 is cer- 
tainly transverse in the generalized Coulomb gauge. The 
right-hand side is also transverse. A mathematical rea- 
son is that the wave equation is found by functional dif- 



Pg(r,i) = 
Mg(r,i) = 



These polarization and magnetization fields are simply 
the sums of the fields Pgm and Mgm produced by the in- 
dividual guest atoms. Finite-order multipole expansions 



ferentiation with respect to generalized transverse func- 
tions. In section IIV Bl it was shown that these deriva- 
tives are always transverse. Physically, the source term 
must be transverse because of current conservation. [Use 
Eq. ^]. 



V. THE QUANTUM MULTIPOLAR 
INTERACTION HAMILTONIAN 

In the previous sections it was shown that the minimal- 
coupling Lagrangian produces the Maxwell-Lorentz equa- 
tions for the electromagnetic fields and the guest charges, 
before and also after choosing the generalized Coulomb 
gauge. Now the goal is to transform the gauge-dependent 
minimal-coupling Lagrangian l|27|l in order to obtain a 
Lagrangian with multipole interaction between the elec- 
tromagnetic field and the guest atoms. The multipolar 
Lagrangian must lead to the same equations for the fields 
and charges, of course. Candidate transformations are 
transformations where a total time derivative of a func- 
tion of the canonical variables is added to the Lagrangian. 
Such transformations leave the action unchanged [l3j . A 
particular transformation of this sort will be used shortly, 
but first some new fields must be introduced. 



A. Polarization, magnetization and displacement 
fields 

In the following, it is useful to describe the guest atoms 
in terms of a polarization density Pg and a magnetiza- 
tion density Mg, rather than in terms of the charge and 
current densities. The former and latter pairs are related 
through [sgf : 

cTg = -V • Pg, Jg = Pg V X Mg. (40) 

In terms of the new variables, the equation of continuity 
is automatically satisfied. The polarization and magne- 
tization fields have the integral representations 



(41a) 
(41b) 



of the polarization and magnetization fields can be found 
by truncating the Taylor expansion in u of the integrands 
on the right-hand sides of the above equations. Such ap- 



/ du {vjnj - Rm) S{r - R„ - u(r„j - R^) ), 
■ Jo 

y2 du w (r,„j - R,„) X f„ij (5(r - R„, - u{r„-,j - R„i) ). 
Jo 
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proximations will be made in section IVDI 

The displacement field D and the magnetic field vector 
H are given by the constitutive relations 

D = eoe(r)E + Pg, H ee /ig - Mg. (42) 

The electric field E and magnetic field B are again de- 
fined by Eq. ||2Jl in terms of a scalar potential $ and a 
vector potential A. It was assumed as before that the di- 
electric is nonmagnetic so that the magnetic permeability 
equals the value /iq of free space. Note that the dielec- 
tric function e(r) is a property of the dielectric alone, 
independent of the guest atoms. On the other hand, the 
displacement field (|42|l does depend on the guest atoms 
because it includes the polarization field produced by 
them. The displacement field defined here is therefore 
different (although the same symbol is used) from the 
displacement field that was defined in section Hi Al where 
no guest atoms were present. In that case, the transver- 
sality of D was evident. The displacement field is also 
transverse when guest atoms are present, according to 
Eqs. (|20c|) and (^UJ in combination with 



B. Classical multipolar Lagrangian and 
Hamiltonian 

The minimal-coupling Lagrangian (|27|l will now be 
transformed to a multipolar Lagrangian by adding to it 
the total time derivative 

-^|drPg(r,t).A(r,i). (43) 

This is the Power-Zienau-WooUey (PZW) transforma- 
tion, and its effect is well-known for free space The 
reason to choose this transformation will be given a pos- 
teriori, when discussing the multipolar Hamiltonian. The 
difference with the free-space PZW transformation is now 
that the vector potential satisfies the generalized rather 
than the usual Coulomb gauge. The transformation was 
already applied to inhomogeneous dielectrics before, in 
[TtI. There it was stated that the polarization density 
Pg in the Lagrangian density should be replaced by a 
"reduced polarization density" at this point, in order to 
stick to the generalized Coulomb gauge for the vector po- 
tential. However, as will be clear shortly, such replace- 
ments are not necessary if functional derivatives with re- 
spect to the generalized transverse vector potential are 
identified as constrained functional differentiations, and 
if computation rules as presented in section Hvl are used 
accordingly. 

After adding the term H43|l . the new Lagrangian imuiti 
can be rewritten as 

-'^muiti = ^ ^"T-crmj ~Vc + y dr £,„uiti- (44) 

mj 



The new Lagrangian density £muiti has the form 
>Cmuiti = ieoe(r)A2--^(VxA)2 + Mg-VxA-Pg-A, 

Z ZflQ 

(45) 

where the definition of the magnetization density (|41b|) 
was used as well as Gauss's theorem. The derivation is 
identical to the free-space case. 

In order to find a Hamiltonian, first the canonically 
conjugate variables must be determined. By reasoning as 
in section lTV CI the constrained functional differentiation 
of the Lagrangian Ly^^m with respect to A produces the 
following field that is canonically conjugate to the vector 
potential: 

n ^ = eoe{r)A - [Pg - eoe(r) V$] = -D. (46) 

Here, the definition of the electric field ((2Jl and the dis- 
placement field H42|) were used. As in the case without 
guest atoms, the field canonically conjugate to the vector 
potential equals minus the displacement field. The dif- 
ference is that now the displacement field also contains 
the polarization field produced by the guest atoms. The 
canonically conjugate field would have been different if 
the minimal-coupling Lagrangian had been used. 

The result Eq. H46() that the canonically conjugate field 
of the vector potential is the full displacement field of the 
medium including guest atoms, is an important general- 
ization of the free-space result In our formalism, it 
could be found rather easily, by realizing that functional 
derivatives must be redefined after choosing a gauge. In 
\it\ . a canonically conjugate field was identified that was 
stated to be different from the displacement field; in 
the matter was reconsidered and the displacement field 
was found as the canonically conjugate field after all, 
but only because the polarization field Pg in the PZW 
transformation H43|l was replaced by a 'reduced polar- 
ization field' for reasons that remain somewhat unclear. 
The effect of the replacement seems to be that func- 
tional derivatives with respect to the vector potential can 
be calculated as partial derivatives, a computation rule 
that in general is valid only before choosing the gauge. 
In contrast, our PZW transformation (|43|l features the 
usual polarization field of the guest atoms, whether we 
choose to do the transformation before or after fixing the 
gauge. We think that our approach is more transparent 
and wider applicable. 

The canonical momenta p„y corresponding to the co- 
ordinate variables q„y of the guest charges are 

where the field 'Pmj stems from the magnetization density 
(|41b|l and is defined as 

F„y = e / Auu B[R,„ - u (r„j - Urn)] x (r„y - R„i). 
Jo 

(48) 
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Note that unUke A and qmj, their canonically conjugate 
variables II and p^j are not fully electromagnetic or fully 
atomic in nature, respectively. 

All canonical momenta have now been determined, so 
that the multipolar Hamiltonian i/muiti can be given in 
terms of the canonical variables (rmj , Pmj) and canonical 
fields (A,n): 

i^multi = ^ Pmj ■ fmj + J dr 11 • A - imulti (49) 



T/(™) I T/(™) 



(50) 



The total Hamiltonian consist of a radiative and an 
atomic part, plus electric and magnetic interactions be- 
tween field and matter. The radiative part of the Hamil- 
tonian is 



rad 



dr 



(51) 



consisting of an electric and magnetic field-energy term, 
respectively. The form of the radiative Hamiltonian has 
not changed after adding the guest atoms, but there is a 
slight shift in its interpretation, since the polarization of 
the atoms is included in the conjugate field. The atomic 
polarization field of atom m interacts with the electro- 
magnetic field as described by 



(m) 



dr 



n 



eoe(r) 



(52) 



In most cases this is the dominant interaction between 
field and matter. The usually weaker magnetic interac- 
tion consists of two terms: 



mj 



(53) 



The first term is linear in the magnetic field and repre- 
sents the paramagnetic energy. Instead of the magne- 
tization Mgm, a reduced magnetization has been 
used in this first term. The reduced magnetization is 
defined as the magnetization [see equation H41bl) ] with 
the Vjnj replaced by Pmj/me The difference has 
been corrected for by a sign change of the second term, 
which is quadratic in the magnetic field. This is the dia- 
magnetic energy of the guest atoms in the nonmagnetic 
dielectric. It can be safely ignored from now on since it 
is much smaller than the other two interactions (see [ll| , 
Sec. 8.6). 

Only the atomic part of the Hamiltonian (|5U|I must 
still be discussed. It has the form 



EPmj 
2m^ 

mj 



dr 



2eoe(r) 



(54) 



The first term in the atomic Hamiltonian H54|l represents 
the kinetic energy of the guest charges; the second term 



is the potential energy of the guest atoms, expressed as 
a polarization energy. The Coulomb term Vc is absent 
in the Hamiltonian, because it cancels against the other 
term quadratic in V<I> that one gets when solving Eq. 146|) 
for A and substituting the result in the Hamiltonian H50|l . 

It is natural to split the polarization energy in Eq. l(5n|) 
into an intra-atomic and an interatomic polarization en- 
ergy, respectively [l^ll^: 



dr 



2eo£(r) 



p2 

2eoe r ^ 

m^n 



P P 

dr 8" ^ 



eo£(r) 
. 

The intra-atomic polarization energy is the potential en- 
ergy that keeps an atom together; the interatomic po- 
larization energy is the only interaction term between 
neutral atoms in the multipolar Hamiltonian l|50() . The 
Hamiltonian is still classical, so that the polarization 
field Pgm (r) (|41a|l is identically zero outside the smallest 
sphere surrounding all charges that make up the (neutral) 
atom TO. For that reason, the classical interatomic polar- 
ization energy (also known as contact energy) is identi- 
cally zero unless bounding spheres of distinct guest atoms 
overlap. In the quantum mechanical description that will 
be given shortly, the expectation value of the polariza- 
tion energy will not be identically zero. The atomic wave 
functions of distinct guest atoms have a non-vanishing 
overlap. However, the overlap falls off exponentially with 
interatomic distance and is negligible unless the distance 
is of the order of the size of the atoms. 

In the rest of this paper, the guest atoms are assumed 
more than a few nanometers apart and their contact en- 
ergies are neglected. Then the atomic Hamiltonian for 
all guest atoms is simply the sum of single-atom Hamil- 
tonians 



o-(m) _ \ ^ "Pmj 



dr 



p2 

2eoe(r) 



The total Hamiltonian Eq. H50|l becomes 



rad 



H. 



(m) 



(m) 



(m) 



M 



(56) 



(57) 



Just like in the free-space case !l3|, in the multipole 
Hamiltonian (|57|l for inhomogcneous dielectrics there 
is no instantaneous interaction term left between well- 
separated neutral guest atoms. This means that in the 
multipolar representation atoms only notice each other 
because they interact with the same (retarded) electro- 
magnetic fields D/[eoe(r)] and B. Of course, the mul- 
tipolar and the minimal-coupling representation should 
give identical physical predictions; in 12,13] the equiva- 
lence is proved for several observables in free space. The 
absence of direct interatomic interactions often makes 
calculations simpler in the multipolar representation. 
This justifies the choice of the PZW transformation H43|) 
out of many candidate transformations. 
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C. Quantum multipolar Hamiltonian 

The goal is now to rewrite the Hamihonian (|57() into 
a second-quantization form, where macroscopic quanti- 
zation has been appUed to the electromagnetic field and 
microscopic quantization to the guest atoms. 

First start with the atomic Hamiltonian, Eq. H5()|l . Fol- 
lowing (standard) quantum mechanics, the electron coor- 
dinates TCmjit) and their canonical momenta Pmj{t) (07)), 
as well as the polarization field Pg, become operators 
that work on the atomic wave functions. The single- 
atom wave functions can be expanded in terms of eigen- 
functions (labelled k) of the single-atom Hamiltonian: 

*m(ri,r2, . . ■,rz^;t) =^Cjnk{t)'ipmk{ri,r2, ■ ■ ■,rz,J- 

k 

(58) 

Second-quantization notation can now be introduced 
by promoting the probability amplitudes Cmk{t) and 
^mki^) become annihilation and creation operators 



with standard anti-commutation relations. These opera- 
tors become the atomic canonical variables in the second- 
quantization picture. (One could even go back and start 
with a Lagran gian that identifies ^'„i and as canoni- 
cal conjugates |22|.') The atomic operators can be written 
as sums over matrix elements. For example, the atomic 
Hamiltonian of atom m in standard second-quantization 
notation is iJ^;"' = Xifc ^^i™^ cl^kCmk- 

The vector potential was again chosen generalized 
transverse and its canonically conjugate field H4ti|) turned 
out to be transverse again, so that these field operators 
can be expanded in terms of generalized transverse modes 
as in Eqs. (|19a|l and l|19bp . The creation and annihila- 
tion operators are written as a\ and a\ now that atoms 
are present. The radiative part H51|) of the Hamiltonian 

becomes i/i-ad = Ea f^xio-l^x + 5). 

The interaction terms Vp"^^ and in second- 

quantization notation become 



(m) _ 



(59a) 



(m) 
M 



J drMg„(r).B(r) = -EE\/^^{«^^rV™fe' / dr M;„^,,, (r) • [V x f,(r)] + H.c.j. (59b) 

A k.k' ^ 



The quantity Pgm.kk' (r) in Eq. H59a|l is the matrix el- 
ement of the polarization field Pg,„(r) with respect to 
states ipmk and ipmk' of atom to. The polarization field 
couples to the field — D/[eo£(r)], which is unequal to 
— E/eo- The interpretation of the interaction is subtle, 
since in the definition (|42|l of the displacement field the 
polarization of the guest atoms is included. The interac- 
tion Ip™'' therefore includes a self-interaction of the po- 
larization field. However, in the expansion in H59a() of the 
interaction in terms of the optical modes, the dielectric 
function £(r) drops out and the coupling becomes rather 
simple. There are no analogous self-interactions in the 
magnetic interaction The magnetic field and the 

reduced magnetization field are canonically independent 
and Mg is not included in the definition of B. All four 
terms in the quantum multipolar Hamiltonian (|57|l have 
now been given in second-quantization notation. 



D. Dipole approximation 

An atom is much smaller than an optical wavelength; 
its spatial structure can not be probed with light. One 
can make the well-known assumption that the polariza- 
tion and magnetization fields associated with the atom 



are concentrated in its center of mass Rm (the nucleus, 
say). Mathematically, this means that the integrands in 
Eqs. (|41al) and (|41bp are approximated by their values in 
u — 0. These two values are the first terms of two infinite 
Taylor expansions in terms of the variable u. The dipole 
approximation is made by keeping only the first term. 

Incidentally, the next terms in the Taylor expansions 
would describe quadrupole interactions, which can be im- 
portant when the guest atoms are not real atoms but 
other (larger) quantum systems in interaction with the 
electromagnetic field. For example, quantum dots ("ar- 
tificial atoms") are much larger than real atoms and so 
their dipole moments can be much larger as well [i^ . 
Quadrupole moments are more important for quantum 
dots than for real atoms, especially when excited in their 
near field by a scanning near-field optical microscope ^3 ■ 
In the following, however, quadrupole and higher-order 
moments are neglected. 

In the dipole approximation, the magnetization (re- 
duced or not) becomes identically zero and the polariza- 
tion field becomes 

Pgm(r) = 5(r - R,„) E cLfe/^il- (60) 

k,k' 

where the atomic dipole matrix elements M^™'' of the 
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^uest atom m are defined as 



(m) ^ 
/^fcfc' — 



-e {ipmk\ y^(r. 



j 



Rm)|V'>i 



(61) 



With equation (|59a|l . it follows that in the dipole approx- 
imation the interaction energy of an atom with the elec- 
tromagnetic field in an inhomogeneous dielectric equals 



(m) 



dip 



E 

kk' 



„(™) 
ik'H-kk' 



D(R™)/[eoe(R 7J 



A kk' 



r 



2eo 



a\c. 



mk 



(rn) 

Cmk' fJ'kk' 



fA(R-n 



H.c. 



(62) 



This gives the important result that inside an inhomoge- 
neous dielectric, a dipole couples to the field — D/[eoe(rll- 
This generalization of the free-s pac e dipole-coupling [ij, 
123 . |23| was also found in 0, 0. As we shall see in 
the following subsection IV EL local-field effects can have 
a strong influence on this interaction. Still, local-field 
effects are often neglected in macroscopic quantization 
theories 



re otten negiec 

EllillilEi 



E. In need of a local-field model 

In section IV CI it was not stressed that the atomic 
Hamiltonian H56|l in general is different for an atom in 
a dielectric and in free space. The potential energy in 
(|56|l . which includes the Coulomb potential, is reduced 
by a factor e(r) as compared to free space. (Such a reduc- 
tion factor is well known for dielectric-filled capacitors). 
As a consequence, energy levels and wave functions will 
be different in a dielectric. This point is missed if one 
starts with a second-quantized description, for example 
when introducing in the medium a "two- level atom" with 
known transition frequency and dipole moment. 

First suppose that the dielectric function is a macros- 
copically averaged quantity that does not change on 
atomic length scales. Then e(r) must be unchanged by 
introducing a guest atom and the Hamiltonian (|56|l can 
be approximated by 



H. 



(m) 



1 

2mp 



E^" 



+ 



drPi 



(63) 



If the guest atom were a hydrogen atom, then its Bohr 
radius and dipole moments would increase by a factor 
e(Rm) and its energy levels would be reduced by the 
same factor, according to this Hamiltonian. Any visible 
line in free space would then be shifted to the infrared 
in a dielectric. The consequences of approximation H63() 
would be that the dielectric has a huge effect on the 
atom's electronic properties. Now the reduction of the 
Coulomb potential (as well as its screening as a func- 
tion of distance) is a well-studied subject in sohd-state 
physics Sometimes one finds the full reduction [as 

described by the Hamiltonian (|63|l ]. while in other cases 
no reduction is found at all. 



In general, outer electronic states of atoms will be more 
affected by the dielectric than the core electrons. An im- 
portant reason for this is the dispersive interaction of the 
guest atoms with the atoms that make up the medium. 
However, such frequency dispersion in the medium is neg- 
lected in the present formalism. We should therefore not 
have the ambition to find an atomic Hamiltonian that 
leads to correct inner and outer electronic states, includ- 
ing medium effects. A modest model is needed that meets 
the requirement that energy levels taking part in the op- 
tical transitions under study should come out right. Such 
a model might be obtained by assuming that the atom 
sits inside an atom-sized cavity with a relative dielectric 
function £(Rm) that is constant inside the cavity; in gen- 
eral e(R„i) will be different both from the macroscopic 
dielectric function just outside the cavity, and different 
from unity (the free-space value). More ambitious de- 
scriptions of medium effects on atomic Hamiltonians re- 
quire at least that dispersion of the dielectric is taken 
into account, perhaps starting from a microscopic model 
of the dielectric 33. 34, 35]. 

A well-known case where reduction of potential energy 
is important occurs when doping solid silicon with phos- 
phorus to make an n-type semiconductor. The high di- 
electric constant of Si (e = 11.7) reduces the potential 
energy between the outermost electron and the rest of 
the P-atom, so that the electron can enter the conduc- 
tion band relatively easily, leaving a P^-ion As 
said before, modifications other than H63|l of the atomic 
Hamiltonian are possible. An important example of the 
other extreme case, where a reduction of the Coulomb 
interaction is absent, will be given shortly. 

We are interested in atomic lifetime changes and line 
shifts caused by the medium. In general, the medium 
induces changes both in the atomic Hamiltonian l|56|) 
and in the atom-field interactions, as compared to free 
space. Effects of the medium that are caused by changes 
in the atomic Hamiltonian will be called electronic ef- 
fects. Changes in atomic dipole moments are an exam- 
ple of electronic effects. On the other hand, effects due to 
modified interactions Fp™'* and Vp^™"* between field and 
atom will be called photonic effects. It is the photonic 
effects, the changes due to altered properties of the elec- 
tromagnetic field, which are of primary interest here and 
in photonics at large. However, only if the electronic 
changes of the atoms are somehow either absent or ac- 
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counted for, can one study the photonic effects. In this 
respect it is fortunate that hue shifts due to changes in 
the interactions (radiative or Lamb shifts) in a medium 
usually are too small to be observable and electronic line 
shifts dominate. Line shifts can therefore be used to es- 
timate medium-induced changes in the atomic Hamilto- 
nian. Given a line shift, one could assign an effective 
dielectric function e(Rm) for the atomic cavity that pro- 
duces the observed transition frequency when inserted in 
the atomic Hamiltonian Eq. 

The distinction between photonic and electronic ef- 
fects is also very important in the interpretation of ex- 
periments. For example, the recently observed fivefold 
reduction of spontaneous-emission rates inside photonic 
crystals is a photonic effect, since possible changes in 
dipole moments were divided out by choosing a reference 
sample with identical electronic effects ^3 ■ Some earlier 
observations of long lifetimes in photonic crystals must 
be attributed to electronic effects [i^ . 

When studying photonic effects of a medium, the guest 
atoms ideally are electronically the same as in free space, 
in particular with the same eigenfrequencies and transi- 
tion dipole moments. In that ideal case, that we refer to 
as the empty-cavity models we have 



£(Rm) = 1 



Vm. 



(64) 



The atomic Hamiltonian is as in the approximation 
Eq. (|63() . now with £(Rm) equal to 1. In other words, 
guest atoms can only be ideal if the dielectric function is 
locally changed to the free-space value 1. The atom-as-in- 
free-space sits inside an empty cavity inside the dielectric. 
A reduction of the intra-atomic Coulomb interaction is 
completely absent in this empty-cavity model. The for- 
mation of such a cavity is beyond the scope of the present 
macroscopic theory. This would require microscopic the- 
ories of the dielectric, involving the Pauli exclusion prin- 
ciple for electrons of both the dielectric and the guest 
atoms. 

The empty-cavity model captures the observed absence 
of large electronic effects of the dielectric on atomic prop- 
erties of interest, but at the same time the model has 
consequences for photonic properties: the local changes 
in £ will give local changes in the mode functions f^, 
and therefore in the dipole coupling Hfi2|l . Atomic 
spontaneous-emission rates will get local-field correc- 
tions. These predictions can be tested experimentally. 
An important example is the emission rate of an atom 
inside an atomic-sized empty cavity in an otherwise ho- 
mogeneous medium. To be precise, £(R) = £ for R 
not coinciding with any of the Rm- The emission rate 
is [3£/(2£ -I- l)]'^^/eTo, where Fq is the free-space emis- 
sion rate jl^ . The well-known in- medium enhancement 
by a factor y/e is further enhanced by the square of a 
so-called local-field factor. Here, the term between the 
square brackets is the empty-cavity local-field factor. 

In a recent study of refractive-index dependent 
spontaneous emission rates, atoms were embedded in a 
low-index molecular complex so as to electronically sep- 



arate them from the medium. For the interpretation of 
the results, it was important that atomic spectra and 
dipole moments did not change appreciably while vary- 
ing the refractive index. The empty-cavity local-field fac- 
tor was indeed observed 'i^. This result is a justification 
for the macroscopic quantization theory for nondispersive 
dielectrics. 

For inhomogeneous dielectrics, it is in general not easy 
to calculate local-field factors, either in the empty-cavity 
model (|64|) or in other models. The simplest assumption 
in the empty-cavity model is that the position-dependent 
local-field factors will have values 3£b(Rm)/(2£b(Rm) + 
1), where £b(Rm) is the bulk dielectric function around 
atom m. The assumption will probably break down when 
£b(R) varies strongly on the scale of the wavelength of 
light. 

The atomic Hamiltonian can be changed in many ways 
and consequently, empty-cavity factors are not the only 
local-field factors that can be obtained from the present 
macroscopic quantization formalism. One could give up 
the macroscopic quantization as being too phenomcno- 
logical and instead describe the microscopic constituents 
of the dielectric in the vicinity of the guest atom. This 
could lead to other local-field factors, depending on the 
question whether the guest atom sits inside a real cavity 
inside the dielectric (of which the empty cavity H64() is 
a sp ecial case), or not. For homogeneous dielectrics, see 
[33L Isol I51I I and references therein. However, for inho- 
mogeneous dielectrics it will be hard to tie a local mi- 
croscopic approach to the macroscopic description of the 
inhomogeneous medium on a larger scale. 



VI. DIPOLE-COUPLING CONTROVERSY 

Many papers appeared in the nineteen-eighties about 
the equivalence of the minimal-coupling and the multipo- 
lar Hamiltonian in free space, for example "20,12112111 
f24j . The Hamiltonians sometimes lead to different results 
in calculations. Some authors argued that the minimal- 
coupling Hamiltonian was to be preferred, while others 
proposed to refrain from using gauge-dependent equa- 
tions to stop the confusion. In the multipolar picture, 
a controversy arose whether a dipole in free space cou- 
ples to minus the displacement field —fx • D/£o, or to the 
transverse part of the electric field —fj, ■ E""". The first 
answer is correct and the book by Cohen- Tannoudji et 
al. helped to settle the argument 0|. It may be use- 
ful to give two sources of confusion even for an atom in 
free space, and to compare the free-space dipole coupling 
with its in-medium generalization Eq. H()2|l . 

The main source of confusion is related to approxima- 
tions. It was found in section IV Dl that a dipole couples 
to minus the displacement field 



-D(R„0/[eoe(Rm)] 



[£(R„OE(R„0 + Pg(Rm)] 



eo£(R-m) 



(65) 
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where Pg is the polarization field of the guest atom it- 
self. In free space, or when assuming an empty-cavity 
model, a dipole couples to the field — D/ep- Still, one 
can find references stating that in free space or in a di- 
electric [l6| a dipole couples to minus the electric field. 
The origin of this mistaken interpretation can be traced 
back to the normal-mode expansion of the displacement 
field operator for the dielectric with guest atoms 



D(r) = ieoe{r) 



A 



2eo 



(66) 



The expanded form of the displacement operator is al- 
most the same as in a dielectric without guest atoms 
[minus Eq. Ijl9b|l ]. but the two differences will now be 
discussed. The first difference, which also plays a role for 
free space, is that the equations of motion of the creation- 
and annihilation operators in (|66|l have terms involving 
the atomic variables, which the equations of motion of 
their counterparts a^"-* and a^^^^ for the field without 
guests do not have. If one approximates the displacement 
field by replacing all the ax and a\ by a^"^ and af^\ 
respectively, then in free space the displacement field H66() 
is equal to the electric field in the absence of the guest 
atoms. In other words, if the guest atoms are taken into 
account in Maxwell's equations, then one finds a dipole 
coupling to the displacement field, whereas a coupling to 
the electric field is found when guest atoms are left out 
of Maxwell's equations. It depends on the observable un- 
der study whether the difference between the two dipole 
couplings can be neglected or not. 

The second difference between the displacement field 
(|^ (with guests) and its counterpart (no guests) is that 
mode functions will be changed locally when guest atoms 
are present, as discussed previously. This difference does 
not show up in the free-space discussion, of course. In 
a dielectric, the above approximation of replacing the 
creation- and annihilation operators docs not make the 
displacement field equal to the electric field. The replace- 
ment would only have this effect for positions r in the 
medium where e(r) equals 1 and where local-field effects 
can be neglected. 

Apart from the main source of confusion, there is an- 
other reason why the interpretation of the dipole cou- 
pling can be confusing: there are two essentially differ- 
ent procedures to go from a minimal-coupling Hamilto- 
nian to a dipole Hamiltonian. The first procedure is to 
rewrite a minimal-coupling Hamiltonian as a multipolar 
or dipole Hamiltonian by canonical transformations of its 
variables (new variables, same Hamiltonian and states). 
The second procedure is a unitary change of picture (new 
Hamiltonian, new states, same expectation values). Un- 
like canonical transformations, picture changes have no 
classical analogues. Confusion is likely to arise when after 
a canonical change a Hamiltonian has exactly the same 
form as after a picture change. The differences between 
the two procedures are excellently presented in |23|| for 
an atom in free space. Either of the two procedures could 



be chosen for dielectrics as well; in this paper, the PZW 
transformation of the Lagrangian was used instead. 



VII. INHOMOGENEOUS MAGNETIC MEDIA 

The present formalism can be generalized to dielectrics 
with inhomogeneous magnetic properties as well, where 
the magnetic permeability becomes fj,o^{r) rather than 
/io. (So here does not represent the magnitude of a 
dipole fi.) Such generalizations are even more interesting 
now that so-called left-handed materials have become 
the subject of intense scientific discussions, after a pre- 
diction that a perfect lens could be made with them [53| . 
Both e(r) and /i(r) of left-handed materials are negative. 
The index of refraction n(r) is also negative and this 
leads to many peculiar properties. Left-handed materials 
will infiuence spontaneous-emission rates of nearby guest 
atoms in different ways than their right-handed counter- 
parts nil. 

Generalizations to magnetic media were already con- 
sidered in 0, and it is relatively straightforward to 
incorporate position-dependent permeabilities in the for- 
malism of this paper, as we will see now. When /i be- 
comes position-dependent, then the only term that will 
change in the classical multipolar Hamiltonian (|57l) is the 
magnetic field energy J dr B^(r)/[2^oM(i")]- The quan- 
tum mechanical description can again be carried out by 
choosing the generalized Coulomb gauge for the vector 
potential. Without guest atoms, the vector potential sat- 
isfies the source-free wave equation 



V X 



1 



■V X A(r) 



0. 



(67) 



The electromagnetic field can be expanded in terms of 
new true modes different from the modes f^. The 
new modes arc the harmonic solutions of the wave equa- 
tion H67|l and so they are generalized transverse, just like 
the modes fA of the nonmagnetic medium. Canonical 
fields can be found by calculating constrained functional 
derivatives for the new Lagrangian in the same way as 
presented in section HVAI It is this point that makes the 
generalization to magnetic media relatively simple. In 
second-quantization notation, the electromagnetic field 
energy becomes huj^ {dldi, -f i), where dj, is the cre- 
ation operator of a photon in the mode (r) . The new 
modes can have mode profiles that differ much from any 
of the modes (r) , but otherwise the theoretical descrip- 
tion of the medium is not much different. In the electric 
and magnetic interactions (|59a|l and l|59b|l of the electro- 
magnetic field with guest atoms, the modes can just 
be replaced by the and the operators c"^^ by d^J^ in 
order to take both the electric and magnetic properties of 
the medium into account. As for nonmagnetic media, in 
the dipole approximation the magnetic interaction (j59b|) 
is zero; the electric dipole interaction H59a|) dominates, 
except for optical transitions with zero dipole moments. 
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When dipole moments are nonzero, the main effect of the 
dielectric becoming magnetic comes from the change in 
the mode functions. 

VIII. SUMMARY AND DISCUSSION 

The aim of the paper was to find a Hamiltonian of 
guest atoms in an inhomogeneous dielectric, with a mul- 
tipolar interaction between the atoms and the electro- 
magnetic field. The multipolar Hamiltonian is simpler 
than the minimal-coupling Hamiltonian, because in the 
former all interactions between the atoms are mediated 
by the retarded electromagnetic field. The main results 
are therefore the quantum multipolar interaction Hamil- 
tonian H57|) together with its dipole approximation (|62|l . 
With this Hamiltonian, one can study how an inhomoge- 
neous dielectric environment can change quantum optical 
processes of resonant atoms. 

In order to obtain the central results, first the elec- 
tromagnetic field was quantized in the absence of guest 
atoms. This has certainly been carried out before 
[T5I ITfil ItR Is^ . but after explaining why real optical 
mode functions can be used whenever e(r) is real, the 
quantization becomes simpler than found in [T^ ITtI l36j . 
In particular, only when real mode functions are cho- 
sen are the associated generalized positions and momenta 
Hermitian [see Eq. H16al) ]. 

Guest atoms were introduced into the theory such that 
Maxwell's equations hold with the atomic charge and 
current densities as source terms. A minimal-coupling 
Lagrangian that gives rise to these equations was eas- 
ily written down. However, after choosing a generalized 
Coulomb gauge, it was not directly clear how to obtain all 
Maxwell's equations in this particular gauge. In section 
IIVI it was explained that the usual functional derivative 
must be replaced by a "constrained functional derivative" 
after choosing a gauge. This is for mathematical reasons 
rather than a matter of taste or convenience. Simple rules 
are given to actually compute these constrained func- 
tional derivatives. As a result, the gauge- independent 
Maxwell's equations were found indeed to hold in the 
generalized Coulomb gauge as well [see Eq. (|39|l ]. 

There is a second advantage of our careful treatment of 
functional differentiation. In the multipolar formalism, 
the field in the dielectric canonically conjugate to the 
vector potential could relatively easily and unambigu- 
ously be identified as minus the full displacement field 
[see Eq. (|46|) ]. This field contains the polarization fields 
produced by both the dielectric and the guest atoms. 

Another important result is that the macroscopic de- 
scription of the dielectric could only be tied up to the 
microscopic description of the atoms by assuming that 
the dielectric function is locally modified by the presence 



of the guest atoms. In particular, in the empty-cavity 
model [Eq. the dielectric function has the value 

1 where wave functions of the guest atoms are nonzero. 
The local modification of the dielectric function will also 
change the dipole coupling, giving rise to local-field ef- 
fects in spontaneous-emission rates. If one would start 
with a two- or three-level description for the guest atoms, 
then one implicitly already assumes a local-field model 
for the dielectric function. Consistency requires to also 
choose that local-field model when calculating the mode 
functions in the dipole interaction (|62|l . 

The application of the Power-Zienau-WooUey transfor- 
mation to the minimal-coupling Lagrangian was shown to 
produce the multipolar Lagrangian, after the generalized 
Coulomb gauge had been chosen. Actually, the gauge was 
chosen earlier than strictly necessary: the choice could 
have been postponed until the canonical momenta were 
determined from the multipolar Lagrangian. The story 
would have been simpler up to that point. The difficulty 
to find all Maxwell's equations would then show up only 
after obtaining the multipolar Lagrangian. The reason 
to first choose the gauge and then do the PZW transfor- 
mation, is that it more clearly shows that the difficulty to 
find all Maxwell's equations was a consequence of choos- 
ing the gauge, rather than a consequence of the PZW 
transformation. 

More generally, the presentation given here is one 
among many possibilities. When going from a classical 
minimal-coupling Lagrangian to a quantum mechanical 
multipolar Hamiltonian, one has to make four steps: one 
step is to choose a gauge, another step is to transform the 
theory to the multipolar formalism. Yet another step is 
made when going from a Lagrangian to a Hamiltonian; 
quantization and second quantization together are step 
number four. These are now given in the order in which 
they occurred in this paper, but the steps can be inter- 
changed. Not all of the 24 permutations are convenient, 
but all routes should lead to equivalent final results. It 
was shown in detail in that step two and three can be 
interchanged for free space: the PZW transformation of 
the minimal-coupling Lagrangian is equivalent to a pic- 
ture change of the minimal-coupling Hamiltonian. The 
equivalence will also hold for inhomogeneous dielectrics. 
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